"Scars" in parametrically excited surface waves 
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FIG. 1. Surface wave patterns in a stadium shaped con- 
tainer at various external frequencies. Adapted from Ref. 7 
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We consider the Faraday surface waves of a fluid in a container with a non-integrable boundary 
shape. We show that, at sufficiently low frequencies, the wave patterns are "scars" selected by 
the instability of the corresponding periodic orbits, the dissipation at the container side walls, and 
interaction effects which reflect the nonlinear nature of the Faraday waves. The results explain 
the observation of a limited number of scars with anomalous strengths in recent experiments by 
Kudrolli, Abraham and Gollub. 



A prominent feature of wave chaos is the "scar" 
phenomenon^. Namely, wave functions are concentrated 
along the short periodic orbits of the underlying classical 
system. In general, the phenomenon appears whenever 
the Eikonal approximation holds, and the corresponding 
rays exhibit ballistic chaotic, dynamics. Since their dis- 
covery by numerical studiesEJ, scars Jaave been observed 
experimentally in microwave cavitieso, tunneling diodesB, 
and capillary wavesu. 

However, theoretical studies of scars! have not ad- 
dressed several factors which may be present in real sys- 
tems, such as dissipation, external forcing, nonequilib- 
rium effects, and interactions associated with nonlinear 
terms of the wave equations. In this Letter we study a 
system where these factors are significant: the surface 
wavesa in a container of non-integrable shape. 

This work was motivated by recent experiments of Ku- 
drolli et al. who_ studied Faraday waves in a stadium 
shaped contained!! The surface waves were excited para- 
metrically by oscillating the container in direction per- 
pendicular to the fluid surface at rest. The experimental 
findings are quite intriguing, see Eig. 1: I. Out of all 
known scars of the stadium billiardEJ only three scarring 
patterns were identified; II. These patterns appear in, 
approximately, 90% of the cases, and their magnitudes, 
compared to the random background, is unexpectedly 
large; III. While probability densities associated with 
eigenstates of the stadium billiard were always symmet- 
ric, nonsymmetric wave patterns have been also observed. 

The purpose of this Letter is to explain this behavior. 
In particular, we show that the excitation threshold for 
observation of a scarring pattern is governed by the Lya- 
punov exponent of the corresponding periodic orbit, bulk 
dissipation, and dissipation at the boundary of the ves- 
sel. These factors suppress most of the scarring patterns 
as well as the random background. The amplitudes of 
those scars which still can be excited are determined by 
nonlinear effects. These effects (interactions) also break 
the symmetry between degenerate scars, which results in 
nonsymmetric wave patterns. 

To begin with, consider the surface waves of an infinite 
system with an infinite fluid depth. Within the linear ap- 
proximation, the amplitude of a plain wave, ake ik r , 



satisfies the equatior 
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where k is the wave number, 7^ = uk 2 is the bulk dissipa- 
tion rate, and v is the viscosity. The angular frequency, 
LUk, satisfies the dispersion relation: 



^t = gk+-k 3 , 
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where g is the acceleration of gravity, £ is the surface 
tension, and p is the fluid density. 

In a finite system, the boundary causes two effects. 
First, it leads to quantization of the wave vector k — ► k a 
(thus also LOk — > w„). When the container is rim full, the 
fluid surface is believed to be pinned to boundarycl. Thus 
the eigenmodes of the surface waves, tp a (r), are solutions 
of the Helmholtz equation in a domain of the container 
shape with Dirichlet boundary conditions. 

The second effect of the boundary is an additional dis- 
sipation due to the viscous flow in the vicinity of the 
container side walls. This dissipation, in a system with 
size L, is of the order y/Uuj/L. For large containers it is 
negligible compared to the bulk dissipation rate. How- 
ever, in small systems and for specific wave patterns, this 
dissipation becomes comparable to 7^. 

Parametric excitation of surface waves, due to verti- 
cal shaking, is accounted by introducing time dependent 
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acceleration, g — > g + Acos(2ujt), into Eq. (|]). The pa- 
rameter space (u>, A) of the resulting Mathieu equation is 
characterized by tongues of instability where ak(i) grows 
exponentially in time. This growth saturates due to non- 
linear effects. 

To take these interactions into account, it is con- 
venient to employ a formalism analogous to second 
quantizations. For simplicity, we first consider an infi- 
nite system. Let us introduce the amplitude variables, 
a k and ak, analogous to boson creation and annihilation 
operators at momentum state k. These variables are o 
numbcrs which can be viewed as boson operators in the 
limit of a large number of particles. The equations of 
motion, then, take the form: 
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where H is the Hamiltonian of the system which consists 
of three terms: H = Ho + Hp + H- m ^. The first compo- 
nent, Hq = X)k w fc a k a k, describes the unforced surface 
waves in the linear approximation, i.e. Eq. (Q). 

The second component, Hp, is the pumping Hamil- 
tonian responsible for the parametric excitation of the 
Faraday waves. Confining our attention to excitations 
near the threshold of the first instability tongue of the 
Mathieu equation, we can write Hp as 



A2ut 



+ C.c. 



(4) 



where A and 2u are the amplitude and the frequency 
of the vibrations. In the quadratic approximation, H ~ 
Ho +Hp, the equations of motion (g) yield 



akjflJJ a e rt , where a 
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+ V/* 2 - (w - ^k) 2 - (5) 



Thus the threshold for instability, at resonance (u> = tuk), 
is h > 7t, . Under this condition, the rate of energy 
pumped into the system exceeds the dissipation rate. 

The third component of the Hamiltonian, H[ n + } is an 
interaction teem. In the mean field approximation it 
takes the forma 

Wint =^4 1) Wkk| 2 |ak'| 2 + 4 2) W4'a-k^ k a_ k , (6) 

kk' 

where I { k j) (9) = T ( k j) (6) + iT k j \d) are complex matrix 
elements which depend on the length of the vectors, 
k = |k| = |k'| as well as the angle, 9, between them. The 
exact fpjm of these matrix elements has been derived by 
Milnert!]. Here we only note that the interaction is local 
and T is of order ujk 2 , while T cx vk A . 

Several mechanisms limit the exponential growth of 
the wave amplitude. They are related to different com- 
ponents of the interaction matrix elements. For instance, 
T k (6) are associated with nonlinear damj DingR while 

T k 2 ^ (9) governs the missphasing mechanismtll. Missphas- 
ing is the phase difference between the pump and the 



excited wave, which appears due to the nonlinearity and 
suppresses the energy absorption. 

The actual limiting mechanism is determined by the 
system size and the excitation frequency. In large sys- 
tems, boundary effects can be neglected, and patterns 
are selected by the local properties of the wave equation. 
This is the regime of pattern formation, where the real 
parts of the interaction terms merely renormalize h and 
Lu k . The dominant—Limiting mechanism, in this case, is 
nonlinear dampinglj. In small systems, contrarily, the 
wave patterns are globally determined by the boundary, 
and missphasing turns out to be more effective. 

A system can be considered as small provided that the 
dissipation time, 1/vk 2 , exceeds the typical time it takes 
a surface wave to cross the system. For capillary waves 
(k > y/ ' gp/Yi) this requirement yields the inequality 



LO < UJ r — 
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In what follows we assume that this condition is satisfied, 
the wave patterns are determined by the boundaxy, and 
missphasing is the dominant limiting mechanisnxj. 

Consider now the Faraday waves in a container of non- 
integrable boundary shape having no discrete symme- 
tries. The modes of the corresponding linear wave equa- 
tion, ^ a (r), are eigenstates of a billiard system of the 
shape of the container. Thus, ipa( r ) are approximately 
random Gaussian functions scarred by the periodic or- 
bits of the underlying classical dynamicsEJ. In principle, 
one can use ip a (r) as a basis for the many body approach 
described above. However, this basis is inconvenient. In 
the experiment, several modes are excited simultaneously 
(approximately 10), and the resulting wave patterns are 
determined by properties of the short time dynamics of 
the system. Therefore, it is more natural to use scars as 
the skeleton of the excited wave patterns. 

Below we develop a "scar representation" . Instead of 
dealing with the amplitudes, n a = |a Q | 2 , of the eigen- 
functions, tp a (r), we consider the amplitudes of scars. 
Using Wigner representation, one can present Ho as 

Ho = )]ii)afia = / e de / dxIU(x;e) n(x) (8) 

a J J 

where x = (r,k) is a phase space coordinate, n(x) is the 
wave amplitude at x, and W(k; e) is the Wigner repre- 
sentation of the imaginary part of the Green functional: 

W(x;e) = ^2S(e-u, a ) J dq e * k -<Va(r-|)C(r+|). 

Assuming the semiclassical limit kL ^> 1 as well as the 
capillary regime, k > (gp/E) 1 / 2 , we can approximate 
W(x; e) as 



W(x;e) ~ <J[e - ff (x)] 1 



]TlU p (x; 



(9) 
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where -ffo(x) = (Sfc 3 //?) 1 / 2 . The first term in the paren- 
thesis is the average of the random excitation over the 
energy shell. The second term is a sum over the primi- 
tive periodic orbits of the system, p. Each term of this 
sum represents the corresponding scar contribution: 



Wp(x;e) ~ 4?rRe 



1 — e --!T+ i (' S W>p) 



*(x p ). 



(10) 



Here S p (uS) = Z p (pw 2 /E) 1 / 3 is the action of the periodic 
orbit with length l p , instability exponent u p , and Maslov 
phase tp p . X p denotes a coordinate on the Poincare sec- 
tion of the periodic orbit. The function S(K p ) is local- 
ized along the p-th periodic orbit with the width of order 
y/L/k. Its integral is unity so that <5(X p ) becomes a 
(^-function in the classical limhxl (i.e. at k — ► oo). 
Substituting (0) in (§), we obtain 



H = fde ed{e)n{e) + ^ u p m) n p , 



(ii) 



where d(e) = ^l(p 2 e/S 2 ) 3 /3 7r is the mean density of 
states in a container of area A, 



= IJ d*S[e-JIo(x)]n(x) 



(12) 



is the average amplitude of the random background, and 



1 



dt n(x p (i)) 



(13) 



is the amplitude of the scar p. Here t p is the period of the 
corresponding orbit, and x p (£) is its coordinates param- 
eterized by the time. The eigenfrequencies of the scars, 
, ,M 



are the poles of (|10|), i.e. the solutions of the equa- 
tion S p {Lo p m ^) = 2irm + tp p — iu p /2. Expanding to leading 
order in u„ we have 
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where A p = u p /t p is the Lyapunov exponent of the orbit 



Now one can substitute the Hamiltonian ( |1 1| ) into the 
equations of motion (^|) and see that the Lyapunov ex- 
ponent plays the role of a dissipation. Indeed, A p /2 is 
the rate at which the "particle escapes from the scar". 
This rate should be added to the dissipation rate. The 
Lyapunov exponents of few short periodic orbits of the 
stadium billiard are tabulated in Table I. Comparison of 
these values with the experimental bulk dissipation rate, 
7b ~ 2 sec -1 , shows that their effect is significant. 

Next, we consider the dissipation at the side walls of 
the container!! 3 ]. This dissipation has been calculated 
for particular cases, e.g., the square and the circular 
containers^! From these results it is evident that bound- 
ary dissipation depends on the particular form of the 
wave pattern. To obtain the result for a general container 



shape, and arbitrary scarring pattern, we first note that 
boundary dissipation takes place very close to the bound- 
ary, at distance of order Id = y/2v/uj ~ 10 -2 cm. View- 
ing the scar as a plain wave scattered from the boundary, 
each collision point adds its own contribution to the dis- 
sipation rate. Thesfi-,contributions can be calculated by 
traditional methodsOJ. Assuming pinning of the fluid sur- 
face at the boundary, the resulting dissipation rate of a 
scar can be presented as a sum over the reflection points 
of the orbit from the boundary: 
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(15) 



where 9 n is the angle between the orbit and the normal 
to the boundary. Examples of 7 P for various scars of the 
stadium billiard are presented in Table I. Note that the 
anomalously large dissipation rate of whispering gallery 
modes (patterns Nos. 8, 9) results from 9 n being close to 
tt/2, see @. 

The boundary dissipation of the random background, 
7, can be obtained from the result for a long ergodic or- 
bit. However, the result diverges since ergodicity implies 
that 9 n can become arbitrarily close to tt/2. Cutting 
this divergence by the wave resolution, A9 n w 1/kL, we 
obtain 

^^l HkL) . (i 6 ) 




X p/2 





3.9 
3.8 
4.8 
3.9 
6.1 
5.9 
6.1 
8.3 



0.35 
0.2 
0.45 
0.4 
0.86 
0.8 
0.68 
2.2 
1.63 



No pattern 




(ffl) 



27 

28.3 

28.7 

29.2 

37.6 

38.1 

49.1 

49.3 

oo 



>/2 



3.4 
4.2 
5.0 
6.2 
2.4 
4.1 
2.4 
3.5 
4.9 



0.42 
0.39 
0.39 
0.8 
0.38 
0.38 
0.36 
0.37 
3.0 



TABLE I. Parameters of some periodic orbits of the sta- 
dium billiard, with dimensions of the experimental system: 
semicircles radius 3 cm, and distance between semicircles 4.5 
cm. l p is the length of the orbit in centimeters, A p is the Lya- 
punov exponent in sec -1 , and 7 P is the contribution to the 
dissipation rate due to the boundary in sec -1 . The velocity 
is set to be similar to the experimental value, 66 cm/sec. 
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Collecting the various dissipation terms and the con- 
tribution of the Lyapunov exponent, the threshold for 
excitation of scars at resonance (a; = Re Wp™') 



is: 



h > h p — 7 b + 7 P + Xp/2 



(17) 



Therefore, large boundary dissipation prevents the ex- 
citation of random Gaussian patterns, as well as other 
modes such as the whispering gallery modes. Moreover, 
increasing the boundary dissipation, e.g. by lowering the 
fluid level, should eventually suppress all of the scars ex- 
cept the horizontal one (No. 2 in Table .0. This is pre- 
cisely what is observed in the experiment! 7 ]. 

Turning to the interaction effects, we first note that 

matrix elements off diagonal in scars, I p L, are usually 

(i) 

much smaller than the diagonal ones, I P p. The reason 
is that scars are localized objects, whereas the interac- 
tion between any two of them is proportional to their 
intersection area. 

Consider, therefore, the scattering from a scar to it- 
self. Here, the wave vectors, k and fc', are either par- 
allel or antiparallel, except for small regions of self in- 
tersections. Thus the diagonal matrix elements of a 
scar associated with an orbit, p, which has a time re- 
versal counterpart — p (such as Nos. 3-9 in Table I) are 



r(2) 



_ 4 2) (0). The 

interaction matrix elements of scars associated with self 
retracing orbits (such as No. 1, 2, and 10 in Table I) are 

= 4f> ~ (4 1} (o) + ifV) + 4 2) (o))/2- 

If only one scar is excited, its amplitude can be evalu- 
ated exactly from the steady state limit of Eq. (||) . For 
scars associated with self retracing orbits, the result is 



where T = Re Ip\p . We assumed that the system was at 
resonance, w = Re lu p , and neglected the imaginary part 

of %L 

Let now two scars be in resonance with the driving fre- 
quency, one can calculate the amplitude of the scar p in 
the presence of the second scar p' by perturbation theory 
in the off diagonal matrix element between these scars. 
The result is 



1 

2T 



h 2 - hi, 



(18) 



,(2) 



,(1) 



(1 - TiQ - T 2 Q 2 



(19) 



where Tj — KeI p J p ]/2T and Q = n^'/rip 1 '. According to 
Eq. (|l^) even a relatively weak interaction can suppress 
the scar with the higher threshold amplitude. Thus one 
should expect to observe only the dominating scar. 

Finally, we comment on the role of discrete symmetries 
in the shape of the container. The exact eigenfunctions 
in this case can be classified according to the irreducible 
representation of the corresponding symmetry group. In 
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the case of the stadium billiard this symmetry is not suf- 
ficient to produce an exact degeneracy of eigenvalues. 
However, the eigenfrequencies of scars, associated with 
an excitation of a set of modes, are degenerate. The 
scars associated with orbits Nos. 3, 4 and 11 in Table 
I can be used as examples. The symmetry between the 
scars is usually weakly broken by perturbations such as 
inhomogeneous pumping or deviation from an exact hor- 
izontal state of the container. Interaction, in turn, drives 
the system into a nonsymmetric wave pattern. 

To summarize, we have shown that Faraday wave pat- 
terns in a container of non-integrable boundary shape 
are mainly scars. The instability threshold for these pat- 
terns is determined by bulk dissipation, Lyapunov expo- 
nents, and dissipation near the container side walls. The 
latter is too strong for random patterns and whispering 
gallery modes to appear. Nonsymmetric wave patterns, 
in a symmetric container, take place due to the existence 
of degenerate scars and interaction effects between them. 
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